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Abstract 

This paper is concerned with a phase field system of Cahn-Hilliard type that is 
related to a tumor growth model and consists of three equations in terms of the vari¬ 
ables order parameter, chemical potential and nutrient concentration. This system 
has been investigated in the recent papers [7] and [9] from the viewpoint of well- 
posedness, long time behavior and asymptotic convergence as two positive viscosity 
coefficients tend to zero at the same time. Here, we continue the analysis performed 
in [9] by showing two independent sets of results as just one of the coefficents tends 
to zero, the other remaining fixed. We prove convergence results, uniqueness of 
solutions to the two resulting limit problems, and suitable error estimates. 
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1 Introduction 

In this paper, we deal with a system of partial differential eqnations related to a model 
for tumor growth that was recently proposed in [16] (cf. also [T7| and [23]) and further 


studied analytically in [71[9l[Tl]. The system reads 

acltyu + - A/i = p((p)((T - 7p) (1.1) 

/r = I3dt(p - A(/? + F'{(p) (1.2) 

dta - Aa =-p{ip){a - ^IJ,) (1.3) 

and it is complemented with the boundary and initial conditions 

dyfi = dyif = duCT = 0 (1.4) 

p(0) = /io, ip{0)=po and a(0) = ao. (1.5) 


Each of the partial differential equations fll.ip - fll.3p is meant to hold in a three-dimensional 
bounded domain endowed with a smooth boundary T and for every positive time, and 
diy in fll.4p stands for the outward normal derivative on T. Moreover, a and f3 are nonneg¬ 
ative parameters, strictly positive in principle, and 7 is a positive constant. Furthermore, 
p is a nonnegative function and F is a nonnegative double-well potential. Finally, po) Uq 
and (Jo are given initial data defined in 

As sketched in [7] (see also its reference list), the physical context is that of a tumor- 
growth model. The unknown function ip is an order parameter which is close to two values 
in the regions of nearly pure phases, say, p ~ 1 in the tumorous phase and ~ — 1 in 
the healthy cell phase; the second unknown p is the related chemical potential, specified 
by fll.2p as in the case of the viscous Cahn-Hilliard or Cahn-Hilliard equation, depending 
on whether > 0 or /? = 0 (see [21[I21[I3]); the third unknown a stands for the nutrient 
concentration, typically, (j ~ 1 in a nutrient-rich extracellular water phase and (j ~ 0 in 
a nutrient-poor extracellular water phase. 

In the paper [5], the asymptotic analysis as both parameters a and (3 tend to zero at 
the same time has been performed, while the cases when each parameter tends to zero 
separately have been left unanswered. The present paper addresses these questions; in 
both cases, the convergence results, the uniquenes for the limit problems, as well as the 
error estimates for the difference of solutions in suitable norms are proved. In particular, 
let us detail here the main mathematical difficulties one encounters in the two passages 
to the limit. 


Passage to the limit as \ 0. The passage to the limit as (3 tends to zero works 
under quite general assumptions on the proliferation function p and on the interaction 
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potential F: p is required to be a nonnegative, bounded and Lipschitz continuous func¬ 
tion, while F is the sum of a convex (possibly multivalued) potential B and a possibly 
nonconvex but smooth part vf. These assumptions are satisfied in quite a large number 
of physically meaningful cases like the classical double-well potential and the logarithmic 
potential, dehned by 

Fci{r) ■■= \{r^ - l)^ = - lyf + \{{l - for r G M (1.6) 

Fiogiy ■= (1 — r) ln(l — r) -f (1 -I- r) ln(l + r) + k{ 1 — r'^y for |r| < 1, (1.7) 

where the decomposition F = B + n is written explicitly. In fll.7p . k is a positive constant 
which does or does not provide a double well depending on its value, and the definition of 
the logarithmic part of Fiog is extended by continuity to ±1 and by -|-cxo outside [—1,1]. 
Moreover, another possible choice is the following: 

F{r) := I{r) + {{1 — r'^y)'^ forrGM, (1.8) 

where I is the indicator function of [—1,1], taking the value 0 in [—1,1] and -|-cxd elsewhere. 
Regarding the function p, in the original model studied in [16] it was set as proportional to 
yF{(p) for |<y5| < 1 and zero elsewhere. Here we can allow such a behavior as well as more 
general ones. Moreover, we can show the uniqueness of the solution to the limit problem 
with a > 0 and (3 = 0, and an error estimate of the expected order 1/2 in (3, under a 
condition of smallness of the hxed coefficient a: the bigger is the Lipschitz constant L of 
the function tt, the smaller has to be a. About this concern, we construct an example of 
severe non-uniqueness for the limit problem fExample 12.41 where aL = 1). 


Passage to the limit as a \ 0. In this second case, when (3 is kept hxed and a tends 
to zero, we obtain similar results, but in a less general setting. Indeed, in this case, due 
to the difhculties in estimating the mean value of the chemical potential p (cf. (14.Ih ). we 
need to assume F to be dehned on the whole of M and to satisfy proper growth conditions 
(cf. fl2.42p L which are, in particular, fulhlled by the classical double-well potential fll. 6 p 
as well as by more general polynomially or exponentially growing potentials. Uniqueness 
and the error estimate can be obtained only under more restrictive conditions on p, which 
must substantially be constant, and F, for which it is required a polynomial growth with 
power four (cf. fll. 6 l) L This is mainly due to the fact that we need to diherentiate in 
time equation fll.2p in order to get these last results (cf. Section 0] for more details). 
However, these results can be compared with [Q] Thms. 2.6 and 2.7] where similar (but 
less restirictive) conditions have to be assumed on F. 


We think that the methods used in our asymptotic analyses could be useful also in 
other situations. In particular, let us point out that in the case of the choice p = 0 
(admitted by our assumptions 02.41) and 02.481) 1 our system 01.1I) - 01.3I) decouples and 
01 . 1 l) - 01 . 2 p reduces to a well-known phase held system of Caginalp type which can be 
seen as a (doubly) viscous approximation of the Cahn-Hilliard system. To this concern, 
let us quote the paper | 20 |, where similar asymptotic analyses are carried out in the case 
p = 0 with the choice 01.61) for F, as both the pameters tend to 0 or just a goes to 0 with 
(3 > 0 hxed. To the best of our knowledge, we do not know of other investigations as 
/3 \ 0: this contribution by us seems to be new also in the case p = 0. Other examples of 
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rigorous asymptotic analyses with respect to parameters intervening on phase held models 
can be found, e.g., in |4H6l[8|[T0 l [TT |[T^ [T9 | [ ^ l2^ . 


Plan of the paper. Our paper is organized as follows. In the next section, we state 
our assumptions and results on the mathematical problem. The last two sections are 
devoted to the corresponding proofs. 


2 Statement of the problem and results 

In this section, we make precise assumptions and state our results. As in the Introduction, 
O C is the domain where the evolution takes place, and T is its boundary. We assume 
n to be open, bounded and connected, and T to be smooth. Moreover, the symbol dy 
denotes the outward normal derivative on T. Given a hnal time T, let 

Q ■= VL X (0,T) and S := T x (0,T). (2.1) 

Moreover, we set for brevity 

V-=H\n), H-=L^{n) and W ■.= {v e ■. dyV = Q onV} (2.2) 

and endow these spaces with their standard norms. For the norm in a generic Banach 
space X (or a power of it), we use the symbol || • ||x with the following exceptions: we 
simply write || • ||p and || • ||* if A = L^iVL) or A = L^iQ) for p G [1, +cxo] and A = V*, 
the dual space of V, respectively. Finally, it is understood that H <zV* as usual, i.e., in 
order that (n, v) = uv for every u ^ H and n G P, where ( • , •) stands for the duality 
pairing between V* and V. 

As far as the structure of the system is concerned, we are given two constants a and 


{3 and three functions p, B and t satisfying the conditions listed below 

a,/5 g(0,1) (2.3) 

p : R —)■ M is nonnegative, bounded and Lipschitz continuous (2.4) 

B ; R —)■ [0, +oo] is convex, proper, lower semicontinuous (2.5) 

T G C'^(R) is nonnegative and tt := vf' is Lipschitz continuous. (2.6) 

We also dehne the potential F : R —)■ [0, +cxo] and the graph i? in R x R by 

F;=.B+7 f and B := dB (2.7) 


and denote by D{B) and D{B) the effective domains of B and B, respectively. It is well 
known that B and the operators (denoted with the same symbol B) induced by B on 
spaces are maximal monotone (see, e.g., [2], Ex. 2.3.4, p. 25]). 

We notice that, among many others, the most important and typical examples of po¬ 
tentials fit our assumptions. Namely, we can take as F the classical double-well potential 
or the logarithmic potential defined in fll.61) and 01.71) . respectively. Another possible 
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choice is the nonsmooth potential fll.8p . In the case of a so irregular potential, its subdif¬ 
ferential is multi-valued and the precise statement of problem fll.lD - fll.Sp has to introduce 
a selection ^ of B{u). 

As far as the initial data of our problem are concerned, we assume that 

fiQ , aQ E H, (fQ E V and ^ , (2-8) 

while the regularity properties we pretend for the solution are the following: 

E H\0,T;V*)nL'^{0,T;V) (2.9) 

(p E H\0,T;H)nL‘^{0,T;W) (2.10) 

^EL\0,T-H). (2.11) 

We notice that fl2.9l) - fl2.10l) imply n,aE C^{[0,T]; H) and p) E C°([0, T]; R). At this 
point, we consider the problem of finding a quadruplet (/i, with the above regularity 

in order that (p, and the related function 


R = - fi) 


( 2 . 12 ) 


satisfy the system 


a{dtp,v) + / dt(pv+ / V/i-Vn= / Rv 
Jo. JQ Jq 

for every v E V, a.e. in (0, T) 

/i = /3dt(p — Aip ^ 7r(y9) and ^ G B{ip) a.e. in Q 


{dta, v) 


Rv 


I VcT ■ Vv = - 
o Jn 

for every v eV, a.e. in (0,T) 

/i(0) = po, V^(0) = (fo and (t( 0) = uo . 


(2.13) 

(2.14) 

(2.15) 

(2.16) 


This is a weak formulation of the boundary value problem fll.ll) - fll.5l) described in the 
Introduction. The homogeneous Neumann boundary condition for p is contained in fl2.10p 
(see 02.21) for the definition of W), while the analogous ones for p and a are meant in a 
generalized sense through the variational equations 02.131) and 02.15p . We notice once and 
for all that the addition of 02.13P and 02.15P yields 


{dt{afi + (f + a),v) + f + a)-Vv = 0 (2.17) 

Jq 

for every v eV, a.e. in (0,T). 

Well-posedness for the above problem is ensured by [9l Thm. 2.2], which states 

Theorem 2.1. Assume 02.41) - 02.7p and 02.8p . Then, for every a,(3 E (0,1), there exists 
a unique quadruplet satisfying 02.91) - 02.lip and solving problem 02.12l) - 02.16p . 


In the same paper [9], the authors study the asymptotic analysis of the above problem 
as both the parameters a and fd tend to zero simultaneously and prove an error estimate 
for the difference of the solution to problem 02.12l) - 02.16l) and the one of the expected 








































6 


Asymptotic analyses for a Cahn-Hilliard type phase field system 


limit problem under further assumptions on the potential F (see [HI Thms. 2.5 and 2.6]). 
Namely, it is assumed that F is everywhere dehned and smooth and satishes some growth 
condition. In particular, the classical double-well potential fll.6p is allowed. 

In the present paper, we discuss the analogous problems obtained by letting just one 
of the parameters tend to zero while keeping the other hxed. The results that deal with 
the possible cases are presented at once. 

In the hrst case, we keep a hxed and let P tend to zero. The limit problem one expects 


is the following: 

{dt{afi + ip),v) + ■ Vv = JqFv VveV, a.e. in (0,T) (2.18) 

=—A(p + ^ + 7r((f) and ^ E B{ip) a.e. in Q (2.19) 

(dta,v) + f^Vcr ■ Vv = — f^Rv VveV, a.e. in (0,T) (2.20) 

(a/j, + (f)(0) = ajuo + (fo and cr(0) = (Tq in , (2-21) 


where R is dehned by fl2.12p . For its solution, we require that 

piEL°^(0,T;R)nL^(0,T;V) (2.22) 

(f E L°°(0,T;V)nL^(0,T;W) (2.23) 

a/i + (f E R^(0,T;V*) (2.24) 

a E R\0,T;V*)nL^(0,T;V) (2.25) 

^EL^(0,T;ff). (2.26) 


Our result on the asymptotic behavior as \ 0 holds under the assumption that a is 
sufficiently small. 

Theorem 2.2. Assume fl2.4p - fl2.7p on the structure and fl2.8p on the initial data. Then, 
there exists oq G (0,1) such that, for a E (0, oq) o^f^d P G (0,1), the unique solution 
(fia,f,Ta,f,o'a,/ 3 ,^a,i 3 ) to problem fl2.12p - fl2.16p . with the regularity fl2.9p - fl2.lip . satisfies 

t^a,i 3 ^ fa weakly star in L°°{0,T] H) n L‘^{0,T]V) (2.27) 

Ta,i 3 ^ Ta weakly star in L°°{0,T]V) r\ L'^{0,T]W) (2.28) 


O'a,y ^ O'a weakly in H^(0,T-,V*) r\ L‘^{0,T]V) (2.29) 

dt(afa,y + Ta,y) -t + Ta) weakly in L^(0, T ; V*) (2.30) 

PTa,y ^ 0 strongly in H^(0,T; H) D L'^(0,T;W) (2.31) 

^a,y ^a weakly in L^{0,T] H) (2.32) 


as P tends to zero, at least for a subsequence, and every limiting quadruplet (fa, fa, o'a, ia) 
solves problem fl2.18p - fl2.2ip . In particular, problem fl2.18p - fl2.2ip has at least a solution 
satisfying fl2.22p - fl2.26p . 

Moreover, we can prove uniqueness for the limit problem and an error estimate under 
an additional restriction on a. 

Theorem 2.3. Assume fl2.4p - fl2.7p and (12.8p . Then, there exists ooo ^ (0, oo) such that, 
for a E (0, ooo), the following conclusions hold: 


i) the solution (fa,Ta,Oa,ia) to problem fl2.18p - fl2.2ip satisfying fl2.22p - fl2.26p is unique; 
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a) if a, 0 )^ 0 , 13 ) a,is the UTiique solution to problem fl2.12p - fl2.16p with the regu¬ 
larity specified by fl2.9p - fl2.lip for fi G (0,1), the estimate 

\\ha,0 — ha\\L‘^{0,T-,H) + ||V^q,/ 3 “ T’a ||l2(0,T;E) + \Wa,0 “ ||L°°(0,T;R)nL2(o,r;E) 

+ ||(«/^o,/3 + y^a,0 + “ (o/ia + + CTo) ||l°°( 0,T;E*) < C'a (2.33) 

holds true for (3 G (0,1) with a constant Ca that depends only on a, R, T, the structure 
of the system, and the norms of the initial data related to assumptions fl2.8p . 

Example 2.4. As said in the Introduction, we can construct an example of severe non¬ 
uniqueness for the limit problem. We take p = 0, so that the third equation is decoupled. 
Moreover, we take the indicator function of [—1,1] as 13 and, given L > 0, we choose tt 
smooth, nonnegative and such that 

|7r'(r)| < L for every r G M and 7r(r) = —Lr for r G [—1,1]. 

Finally, we take 

7^0 = ho = o'o = 0. 

Then, for any function G L°°(0,T) satisfying \'fi{t)\ < 1 for a.a. t G (0,T), the definition 

fi{x,t) :=—Lfiit), ip{x,t) := fi(t), a{x,t) := 0 and ^{x,t) := 0 

provides a solution if aL = 1. Indeed, the initial and boundary conditions are trivially 
satisfied as well as the equations, since we have 

a/i if = —aLfi = D, fx = —Lfi = —Lip = 7 r{p) and ^ = 0 G = B{p). 

In the second case, we keep (3 hxed and study the asymptotics with respect to the 
parameter a and the corresponding expected limit problem, namely 

{dtPiV) + p-Vv = /qHu VveV, a.e. in (0,T) (2.34) 

/u =/3dtp — Ap-h f 3 - 7 r(p) and ^ E B{p) a.e. in Q (2.35) 

(dtcr,v) 3 -JqVct ■ Vv = — f^Rv VveV, a.e. in (0,T) (2.36) 

p(0) = po and (t(0) = cxo in R, (2.37) 

where R is dehned by fl2.12p . For its solution, we require the following regularity: 

fiE L^(0,T;V) (2.38) 

p E Il\0,T;H)nL\0,T;W) (2.39) 

a E H\0,T;V*)nL\0,T;V) (2.40) 

^EL\0,T;H). (2.41) 

Our results hold in a less general setting. We need a hrst restriction on the potential F in 
order to prove an asymptotic result and the existence of a solution to the limit problem. 
Namely, we also assume that 

D(B) = R and \B^{r)\ < CB{B{r) 3- l) for every rGM, (2.42) 

where B^{r) is the element of B[r) having minimum modulus and Cb is a given constant. 
Notice that the classical potential fll.6p and similar potentials with polynomial or expo¬ 
nential growth satisfy this assumption. Moreover, fl2.42l) still allows B to be multi-valued. 
We have the following result: 
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Theorem 2.5. Assume fl2.4p - fl2.7p and 02.421) on the structure and 02.8p on the initial 
data. Moreover, fora, (3 G (0, 1), let the unique to problem fl2.12p - 


02.16P satisfying 02.9p - 02.1ip . Then we have that 

Ta,g h-y weakly in L‘^{0,T-,V) (2.43) 

Ta,0 ^ Ty weakly in H) n L‘^{0,T-,W) (2.44) 

weakly in H^{t),T-,V*) A L'^{t),T-,V) (2-45) 

^ C/3 weakly in L^(0, T; H) (2.46) 

OiP'a,i 3 ^ 0 weakly in H^{0,T;V*) and strongly in Lf if), T-,V) (2.47) 


as a tends to zero, at least for a subsequence, and every limiting quadruplet {fig, ipp, ag, C/?) 
solves problem 02.34p - fl2.37p . In particular, the limit problem fl2.34p - 02.37p has at least a 
solution satisfying 02.38p - 02.411) . 

Both the uniqueness for the limit problem and the error estimates look difficult to 
prove. We can overcome such a difficulty only in a particular case. Moreover, we have to 
take p constant and make stronger conditions on the potential F (which, however, still 
allow the choice of the classical potential Ol.bp i and slightly reinforce those on the initial 
data. Namely, we also assume that 

p is a nonnegative constant (2.48) 

F is a function on R satisfying |F"(r)| < C(r^ + 1) 

for every r G R and for some constant C > 0. (2.49) 

We observe that in this setting the last condition in fl2.8p is a consequence oi ipo and 
fl2.49p . Indeed, these assumptions ensure that F{r) = O(r^) as |r| tends to +cxo and that 
<Po £ thanks to the Sobolev inequality. 

Theorem 2.6. Assume 02.41) - 02.71) . 02.421) and 02.481) - 02.49p on the structure, and 02.8p 
on the initial data, and let (3 G (0,1). Then, the following conclusions hold true: 

i) the solution {fig,(pg,ag,^g) to problem 02.34p - 02.37p satisfying 02.38p - 02.411) is unique; 

ii) if is the unique solution to problem 02.12p - 02.161) with the regu¬ 

larity specified by 02.91) - 02.Ill) for a G (0, 1), the estimate 

\\Ta,y — Ty\\L‘^(o,Ty) + ||7’o,/3 — Ty\\L^{0,T-,H)nL‘^(0,T-y) 

+ \Wa,g — <7/3||L°°(0,r;H)nL2(0,T;Y) — (2.50) 

holds true for a G (0, 1) with a constant Cg that depends only on (3, hi, T, the structure 
of the system, and the norms of the initial data related to assumptions 02.8p . 


Now, we list some facts. We repeatedly make use of the notation 

Qt := II X (0,f) for t G [0,T] 

and of well-known inequalities, namely, the elementary Young inequality 
ab < 5a^ ^ every a, 5 > 0 and 5 > 0 , 


(2.51) 


(2.52) 
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the Holder inequality, and its consequences. Moreover, as R is bounded and smooth, we 
can owe to the Poincare and Sobolev-type inequalities, namely, 

Iblly < C^IIVpIIh + for every e R (2.53) 

V C L^(R) and HuHq < C|ln|ly for every v e V and 1 < g < 6 (2.54) 

L'^(R) C V* and ||n||* < C||n||g for every v G L^(f2) and q > 6/5 . (2.55) 

In fl2.53p - fl2.55p . C only depends on R. Furthermore, we denote by A the Riesz isomor¬ 
phism from V onto V* associated to the standard inner product of V, i.e., A \ V ^ V* is 
dehned by 

{Au,v) ■.= {u,v)y = / (yu-'Vv + uv) ioTU,v^V. (2.56) 

Jn 

We notice that Au = —Au + u ii u E W. We also remark that 

{Au,A~^v*) = {v*,u) for every u eV and v* E V* (2.57) 

{v*,A~^u*) = {u*,A~^v*) = for every u*,v* E V *, (2.58) 


where (•, • )* denotes the dual scalar product in V* associated to the standard one in V, 
and recall that {v*,u) = fQV*u if v* E H. As a consequence of fl2.58p . we have 


d_ 

dt 



2{dtV*,A ^v*) ioY every V* E 


(2.59) 


Finally, in order to simplify the notation, we follow a general rule in performing our a priori 
estimates. The small-case italic c without any subscript stands for different constants, 
which may only depend on R, T, the shape of the nonlinearities and the norms of the 
initial data related to assumptions at hand. A notation like cs signals a constant that 
also depends on the parameter 6. We point out that c and cs do not depend on a and (3 
and that their meaning might change from line to line and even within the same chain of 
inequalities. On the contrary, constants that are later referred to are always denoted by 
different symbols, e.g., by a capital letter. 


The starting point for the proofs given in the next sections is the following general 
result (see P Thm. 2.3]): 

Theorem 2.7. Assume fl2.4p - fl2.7p and fl2.8p . Then, for some constant C which only 
depends on R, T and the shape of the nonlinearities, the following holds true: for every 
a,(3 E (0,1), the solution to problem fl2.12l) - fl2.16l) with the regularity specified 

by fl2.9l) - fl2.1ip satisfies 

+ II V/i||L2(o^T;R) 

+ + \\t\\l<=°{0,T-,V) + l|-^(‘d)llL“(0,r;Ll(O)) 

+ ||<9i(a/i -f ^)\\l^{0,T-,V*) + ||cr||Rl(0,T;y)nL°°(0,T;R)nL2(0,T;y) 

< C (o^/^IIp-oIIr + IIt’oIIv + ll-^(‘^o)|liqQ) + IkollR) (2.60) 

as well as 


||h||L2(0,r;y) + \\t\\l'^(0,T-,W) + ||'Cl|L2(0,r;R) 

< C (a^/^||/io||R + IIt^oIIv + \\F{To)\\l^q) + Ikolln + ||h||L2(o,r;R) + l)- 


(2.61) 
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3 Proofs of Theorems 12.21 and 12.3 


It is understood that the assumptions fl2.4p - fl2.7p and fl2.8p are in force. 

Proof of Theorem 12.21 We start from the uniform estimates stated in Theorem 12.71 
As a is hxed, fl2.60p ensures that even is bounded, so that fl2.6ip provides 

a bound for its left-hand side. Therefore the convergences specified in fl2.27p - fl2.32p hold 
for a subsequence. Moreover, as '■= <yiia,i 3 + ^a,i 3 converges to r]a '■= Oifia + weakly 
in 77^(0, T; l/*)nL^(0,T; V), by applying the Aubin-Lions lemma (see, e.g., [HI Thm. 5.1, 
p. 58]) we deduce that 


Va,i 3 —t rja strongly in L^(0,T;77) as /I \ 0. (3.1) 

Now, we prove that fl3.ip implies that both 

and stTOUglj m L‘^{0,T; H) (3.2) 


provided a is small enough. We show fl3.2p by using a Cauchy sequence argument: we 
write fl2.14p for the solutions corresponding to /5 and f3' and take the difference; then, we 
multiply by a, sum ipa,i 3 — ^a,/ 3 ' to both sides, rearrange and hnally test by ipa,i 3 — 

We obtain 


/ ((ha,/3 - Va,l3') - a(/5 dtipa,l3 “ dtiPa,IS')) 

Jn 

ll7^a,/3 f7Q:,/3'|ljr ~\~Oi |V(9?a,^ ^a,/3')\ 




a.e. in (0,T). Then we integrate over (0,T) and observe that the resulting left-hand 
side tends to zero as (3,(3' \ 0, due to the strong convergence given in fl3.ip and fl2.3ip 
coupled with the weak convergence fl2.28l) . The term a — ^a,/ 3 ')i^a,i 3 — ^a,/ 3 ') is 

nonegative thanks to the monotonicity of 7?; the last term can be treated using the 
Lipschitz continuity of vr, namely. 


Jq 

where L denotes a Lipschitz constant for vr. Hence, from fl3.3l) and the subsequent remarks 
we deduce that is a Cauchy sequence in L‘^{Q) = L^(0,T;77) provided aL < 1, so 

that fl2.28p entails that (pa ,0 converges to strongly in L^(0, T; 77) and fl3.2l) is completely 
proved. 

From fl3.2p it follows that 7 r{(pa, 0 ) and p{(pa, 0 ) converge to 7i{(pa) and p{(Pa), respec¬ 
tively, strongly in L^(0,T;77). Therefore, we can identify the limits of the nonlinear 
terms ^a ,0 and Ra, 0 - For the former, we can apply, e.g., [H Cor. 2.4, p. 41]. For the latter, 
we note that Ra ,0 converges to p{(pa){o'a — Pa) strongly in L^{Q) since (cf. 02.291) ) we 
also have a strong convergence of to Uq in L‘^{Q). At this point, we can write the 
integrated-in-time version of problem fl2.13p - fl2.15p for the approximating solution with 
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time dependent test functions and take the limit as a tends to zero. We obtain the anal¬ 
ogous systems for and this implies fl2.18p - fl2.2ip for such a quadruplet. 

□ 

Proof of Theorem 12.31 In order to show Theorem 12.31 we do not follow the order of 
the statement. Indeed, assume for a while that the part ii) has been proved for every 
solution to problem fl2.18p - fl2.2ip (with a constant Ca that might depend on the solution 
we are considering), provided that a is small. From this, we derive the uniqueness part i). 
Indeed, let i = 1,2, be two solutions. If a is small, inequality fl2.33p holds 

for both of them with a common constant Ca- Hence, we immediately derive, e.g.. 


Ilhl ~ h2||l,2(o,T;R) + ||</^2 “ </52 || L2(0,T;V) 

+ Iki - cT 2 ||L 2 (o,r;V) < 2C„/3^/^ for every j3 e (0,1). 


This implies that pi = p 2 , = ^2 and ai = 02 - Then, = ^2 by comparison in fl2.19p . 

Now, we show that the error estimate fl2.33p holds for every solution (/ia, ipa, o'a, ^a) to 
the problem fl2.18p - fl2.2ip (with a constant Ca that might depend on the solution we are 
considering). To this end, we present equations fl2.14p - fl2.15p . fl2.19p - fl2.20p . fl2.17p and its 
analogue obtained by summing fl2.18p and fl2.20p in a slightly different form. Namely, in 
each equation, we add the same function to both sides and make the Riesz ismorphism A 
appear (see fl2.56p ). In doing so, we owe to the regularity conditions fl2.9p - fl2.lip and 
fl2.38p - fl2.4ip . If {ji,Tp,'a,^ stands for the solution to the limit problem, we have 


dt{apa,/3 + V^o,/3 + fLa,/?) + -^(ho,/? + <^ 0 ,/?) — l>‘a,l3 + 

i^9t(pa,l3 T Aipa,l3 T ^ce,0 T 
9t(^a,l3 T A,a,0 T ^a,/3 

dt{ajl + Tp + a) + Aijl + a) = p -h a 
71 = ATp -h ^ -h vr(^) - Tp 
dfC Ag = —R -|- G, 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 


where Ra,^ and R are dehned by fl2.12p according to the equations we are considering. 
All these equations are meant in the framework of the Hilbert triplet {V, H,V*), a.e. 
in (0,T). However, the explicite versions of fl3.5p and (13.8p also hold a.e. in Q. Moreover, 
we have to add the conditions ^a ,/3 G B{,^a,p) and ^ G Bijp) a.e. in Q as well as the 
initial conditions fl2.16p and fl2.2ip . respectively. Now, we take the differences between 


fl3.4p - fl3.6p and fl3.7p - fl3.9p and have 


dt{ap + p + g) + A{p + g) = p + G 
P = (5dtPa,0 +Ap + ia,0 - ^ + 7r(</?a,/3) - 7r(^) - p 
diG -\- Ag = —{Rap — R) T ^ 1 


where we have set, for convenience, 

/i . Pap pj ■ ^ap and G . ^aP 

At this point, we write these equations at time s G (0,T) and test them by 

A~^{ap + p + g){s), —p{s) and cr(s), 
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respectively. Next, we sum up and integrate over (0,t) with respect to s, for an arbitrary 
t G (0, T). Then, we recall and use fl2.56p - fl2.59p : by rearranging a little and omitting the 
evaluation at s inside integrals for brevity, we obtain that 


\\{afi + ip + a){t)\\l + / {fi + a){afi + ip + a) + 


V ds 


' Qt 


+ / 

jQt 


^ Qt 




= / {fi + a,ap + (p + a)^ds - P / dtipa^pip+ / {(p - (7r((p„,^) - 7r((p))}(p 


^ Qt 


^ Qt 


- / {Ra,0 - R)o- + / |a| 

^ Qt Qt 


(3.10) 


The terms on the left-hand side not having a dehnite sign are treated simultaneously in 
the following way: 


J Qt 
> a 


{p + a){afi + ip + a) - fiip= [a\p\‘^ + pa + a{afi + ip + a)) 


^ Qt 


' Qt 


pa 


> 


a 


' Qt 


' Qt 


' Qt 


\a\ 


\a\\v \\ap + p} + (t||* ds 


\a\\yds — C5 / llap-I-(p-I-( t||^ ds. 


^ Qt 


Now, we deal with the right-hand side of fid.iop . For the first integral, we observe that 
the norm of the embedding id C Id* is 1 (since the norms of V and H are the standard 
ones), and we have the estimate 


a 


(p + a,a/i + (p-l-cr)^ds < - 


+ Ill'll*) ds-I-c / ||a/i-I--I-cr||^ ds 


< 


a 


+ 


a 


\a\'^ + c / ||a/i-I- 9 ?-I-(j||^ ds 


^ Qt 


^ Qt 


Next, we have that 


/S / dtp>a,pp><5 / |(p|^ -1- / \dtiPa,p\‘^ < S 


y ds cs/d , 


' Qt 


^ Qt 


^ Qt 


the last inquality being due to fl2.60p for dt<Pa,p- For the next term, we use the Lipschitz 
continuity of vr and still denote by L the Lipschitz constant of tt. Then, we obtain that 


{lP - {7T{(Pa,f}) - Tr{p))) }ip<{l + L) / 

£ Qt 

(1 + L) / ^{ap + ^ + a) — a (pfi — / Lpa I 
\J Qt J Qt J Qt J 


< 5 


+ cs \\ap + (p + a\\l ds 


a 


' Qt 


\p\^ + {l + Lfa 


ds + 6 


y ds + Cs 


\<y\ 


^ Qt 
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and we assume at once a to be small in order that (1 + L)^a < 1. Finally, there is one 
more term to treat on the right-hand side of fl3.10l) . namely, 

- {Ra,^-R)a=- 

jQt 


' Qt 


< 


' Qt 


-P{P)\ Wa,/3- Pa,/3\ |o-| 


\p{p)\ |a-/i| |a| 


^ Qt 


<c \(p\\(Jot,0 - Pot,f}\W\ +C \a-fi\\a\ 


^ Qt 


^ Qt 


Note that the last inequality holds true since p is Lipschitz continuous and bounded. On 
the other hand, we can use the Holder and Sobolev inequalities to obtain 


^ / Iv^l h'cti/Sl 1^1 — ^ 


4 11^0,/I /^a,/3||4 ||^||2 d,S 


' Qt 


< C 


V lka,/3 — Pa,p\\v ||cl||r < 5 


ds ~\~ I ||cra,/3 /^q:,/ 3|1\/ ll^llir ds 


and we observe at once that the function s IKca,/? ~ ha,/ 3 )('S)|lv is bounded in L^(0,T), 
thanks to fl2.60p for aa,j 3 and fl2.60p - fl2.6ip for pa,i 3 - Finally, we have that 


/ \a-p\\a\<^ f |/ip + 

jQt ® JQt 




^ Qt 


At this point, we collect fl3.10p and all the inequalities we have obtained, rearrange and 
infer that 

\\\{ai, + ^ + a)ml + z f 


-|- (l — (1 + R)^o. — 4(5^ 


' Qt 

ft 


V 


ds -\- 


|cr(t)p + {1 - 6) / ||(T||y ds 


<C5 {l + \\(Ja,p - PaMlv) IWUnds + cs / \\ap + ifi + a\\l ds + Cs(3 . 

Jo Jo 

Then, we choose <5 small enough (let us recall that we are assuming (1 -|- L^a < 1) and 
apply the Gronwall lemma. This yields fl2.33l) . and the proof is complete. □ 


4 Proof of Theorems 12.51 and 12.6 


It is understood that assumptions fl2.4p - fl2.7p and fl2.42l) on the structure and fl2.8l) on 
the initial data are in force. We hrst prove Theorem 12.51 The main tool is Theorem 12.71 
applied to the solution {pa,/ 3 , (Ja,/ 3 , ia,p) to problem fl2.12l) - fl2.16l) . However, we need a 

preliminary estimate that has already been performed in [9]. Nevertheless, for the reader’s 
convenience, we repeat the core of the argument here. 

Auxiliary a priori estimate. We omit the indices a and (J for a while. We observe that 
the growth condition fl2.42p formally implies (see [9l formulas (2.25)-(2.26) and Rem. 2.5]) 

[ < c f [B{Lp{t)) + l) for a.a. t e (0,T). (4.1) 

Jn Jn 
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Now, we simply integrate fl2.14p over and use the homogeneous Neumann boundary 
condition for 9 ?. We deduce that 

[ /i(t) = f (3dt^p{t)+ [ i{t)+ f for a.a. t G (0,T). 

Jn Jn Jn Jn 

As TT is Lipschitz continuous and vf is nonnegative, the above identity and fl4.ip imply that 


/i(i) < / P\dt^{t)\ + c / F{^{t)) + c / \^{t)\+c. 


(4.2) 


By taking advantage of fl2.60p . we deduce that the function t i-A is estimated 

in L^(0,T). Therefore, by using fl2.60p once more and the Poincare inequality fl2.53p . we 
conclude that 

ll/^l|L2(o,r;V) < c- (4.3) 


Proof of Theorem 12.51 Now, we explicitly write the indices a and jS. First, fl4.3p 
ensures that the left-hand sides of both fl2.60p and fl2.6ip are bounded, so that, since f3 is 
hxed, we deduce fl2.43p - fl2.47p . at least for a subsequence: note that fl2.47p follows from 
02.43p . f|2.44p and the bound for in L^(0, T; 1/*). Now, let (yU/ 3 , y?/?, a/?, ^/j) 

be any limiting quadruplet. Then, fl2.44p - fl2.45p imply that the initial conditions for 
(v^/ 3 , P"^) are satisfied. Moreover, the Aubin-Lions lemma (see, e.g., [IHl Thm. 5.1, p. 58]) 
ensures that (fa ,0 converges to (pp strongly in L^(0, T; H) (even better, of course), whence, 
71 ( 7 ^ 0 ,/?) and p{(Pa,/ 3 ) converge to 7r((^^) and p^P/s), respectively, strongly in L^(0,T;i7). 
Therefore, we can identify the limits of the nonlinear terms and Rap- For the former, 
we can apply, e.g., IH Cor. 2.4, p. 41]. For the latter, we note that Rap converges to 
p{pp){aj 3 — Pjs) (at least) weakly in L^{Q). At this point, we can write the integrated-in¬ 
time version of problem fl2.13p - fl2.15p for the approximating solution with time dependent 
test functions and take the limit as a tends to zero. We obtain the analogous systems for 
(yW/?, (T/ 3 , ^^), and this implies fl2.3411 - 02.361) for such a quadruplet. □ 

Now, we assume that fl2.48p - fl2.49p hold, in addition, and start proving Theorem 12.61 
However, as in the previous section, we do not follow the order of the statement. Indeed, 
let us assume for a while that its part ii) has been proved in the following modihed version: 
the error estimate fl2.50p holds for every solution to the limit problem fl2.34l) - fl2.37l) (with 
a constant Cy? that might depend on the solution we are considering). From this, we 
derive the uniqueness stated in i) in the following way: let (/ij, pi, ai, ^i), i = 1, 2, be two 
solutions. Then inequality fl2.50l) holds for both of them with a common constant Cp. 
Hence, we immediately derive, e.g., 

Ilhl ~ h2||L2(o,T;y) + IIT’I — P2\\l^{0,T-,V) + Ho'l “ O'2\\l^(0^T;V) < Cp 

for every a G (0,1). We deduce that pi = /i 2 , pi = P 2 and ai = a 2 , whence also .^1 = ^2 
since B is single-valued. 

So, it remains to prove that the error estimate fl2.50l) holds true for every solution to 
the limit problem fl2.34l) - fl2.37l) (with a constant that might depend on the solution we 
are considering). To this end, we need a further estimate on the solution (pap, Pap, ^ap) 
to problem fl2.12p - fl2.16p . 
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A new a priori estimate. We proceed formally in order not to be too technical. 
Moreover, let us omit the indices a and /9. We consider equation fl2.13p . as well as the 
variational equation one obtains by formally differentiating fl2.14p with respect to time. 
Then, we present them in the form of abstract equations by introducing the Riesz operator 
A defined by fl2.56p . i.e., 

dtiptii + (p) + Afi = p{a - /i) + /i 

dtp = I3dtp> + Adtp> + F"{p>)dtp> - dt^p. 


These equations are meant in the sense of the Hilbert triplet {V, H,V*), a.e. in (0,T). 
Now, we test them by A~^dtp and —A~^dt(p, respectively: we sum up and integrate 
over (0,t). By omitting the evaluation point inside integrals for brevity, we obtain 

a f {dtp, A~^dtp) ds+ f {dt(p,A~^dtp) ds+ f {Ap,A~^dtp) ds 

Jo Jo Jo 

+ p I {p, A~^dtp) ds - I {dtp, A~^dt(p) ds + (3 j {dfip, A~^dtip) ds 
Jo Jo Jo 

+ f {Adtp>,A~^dtp>) ds + [ {F"{ip)dt’p,A~^dt’p)ds— f {dt’p,A~^dt’p)ds 

Jo Jo Jo 

= p {a,A~^dtp)ds+ / {p,A~^dtp) ds. 

Jo Jo 

Now, let us recall fl2.56l) - fl2.59p . Two terms cancel out by (12.581) and we have 


= -/ {F"{(p)dt(p,dt(p)^ds + 


ds + p (a, dtp)^ ds 


5 IIa‘(<)K -1 ll/<oli; +1 ll/ioll?, +1 IIa^oII; + f II9.4 ’(o)II: 


(4,4) 


Therefore, we just have to deal with the right-hand side. In order to control the integral 
involving F", we account for the dual Sobolev inequality fl2.55p . the growth condition 
fl2.49p . the Holder and Sobolev inequalities, the estimate fl2.60p for (p, and have 

- f {F"{ip)dtip,dtip)^ds < f \\F”{ip)dtip\\4dM\*ds 


< c 


\F'\ip)dtip\\6/5\\dtip\\^ds < c / \\1 + ip'^\\3\\dtip\\2\\dt(p\\*ds 


<C (l + llv^lle) WdMUds <c (l + \\ip\\lr) \\dtip\\H \\dt^\\*ds 




'WWh 


ds + c 


I ds. 


Another term to treat is 

p {a,dtp)^ds =-p {dta,p)^ds + p{a{t),p{t))^-p{ao,po)^ 


< 


p 


p 


ds + p\\dta\\L2^o,T-y*) + tII/^WIL + P lkllL-(o,T;y) + C 
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the last inequality being due to fl4.3p (which holds true under the assumptions of Theo¬ 
rem |231 thus also in the present case) and fl2.8p . The next term on the right-hand side 
of fl4.4p is (l/2)||p(f)||^. We observe once more that the norm of the embedding H <zV* 
is 1 since we are using the standard norms in V and H. Therefore, we have 

Finally, we consider the norms of the initial values of the time derivatives. We for¬ 
mally have 

f3dt'^{0) = po — A(po + ~ ) 

which is a hxed element of V* due to 02.81) . Indeed, for F'((Pq), we make the following 
observation. Assumption 02.491) implies \F'{r)\ < c{\r\^ + 1) for every r G M. As (po £ V, 
we have (po ^ by the Sobolev inequality. We infer that G whence 

F'{(po) G H . At this point, we can collect 04.41) and all the inequalities we have proved, 
then apply the Gronwall lemma. We conclude that 

<y^^‘^\\dtPa,/3\\L^(0,T;V*) + P^^‘^\\Pa,/3\\L°°(0,T;V*) + ||1|L°°(0,T;V*) < C, (4.5) 
where we have used the full notation with indices. 


Conclusion of the proof of Theorem 12.61 We prove that 02.501) holds for every 
solution (/i^, (T^) to problem 02.34l) - 02.37l) (with a constant Cjs that might depend on 

the solution we are considering). We often omit writing the evaluation point, explicitly, 
in order to simplify the notation. We take the difference between 02.13l) - 02.15l) . written 
for (/ia,/ 3 , da,^), and ( j 2.34 | ) -(|2^, written for {np, pp, (Tfi), where and 

^j 3 = By setting, for convenience, 

fi . Pa,0 P0 1 P • Pa,0 P0 and (J . ^a,0 ^0 ) 

and recalling that p is a nonnegative constant (cf. fl2.48p h we have 


/ dtpv+ / Wp-Vv = -a{dtiJ,a, 0 ,v) +p / {a - p)v 
Jn Jn Jn 

/ pv = (3 / dt(pv+ / V(p-Vv+ / (Ca,/3 - ^/3 )t- h / {'K{Lpa, 0 ) - Tl{p 0 ))v 
Jq Jfi Jq Jq Jq 

{dta,v) + / Va -Vv = p / {p - a)v , 

Jn Jn 

where each equality holds true for every v ^ V and a.e. in (0,T). Now, we take the test 
function v equal to 

if -|- Pp, p — p and a , 

respectively. Next, we integrate with respect to time, exploit the initial conditions 
(p(0) = 0 and cr(0) = 0, add the volume integrals of P/9|pp and of p\a\^ to both sides 
for convenience, and rearrange a little. We have, for every t G [0,T], 

1 


\p{t)\^ + l3 dtpp+ Vp-Vp + M |V/i|^ 


' Qt 


' Qt 


^ Qt 


(l+P/9) 


/S 


' Qt 


r + |/ / d,pp+ I iVpp 


Vp • V/i 


'O 


+ / {^a ,0 — ^ 0 )P + 


^ Qt 


J Qt Qt 

\a{t)\^+ [ \Va\^+p[ |a| 
^ Qt ^ Qt 


' Qt 
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J Qt J Qt J Qt 

+ |3^J) ds+ p j {{a - p){p + (dp-a) + I3\p\^+ \a\^]. (4.6) 

JQt 


a 


Four terms on the left-hand side cancel out, and the other ones are nonnegative. Now, 
we estimate each integral on the right-hand side, separately. For the hrst two of them, 
we have ^ 

/ ^ j Ihr + c/ \p\^ 

JQt JQt ^ JQt JQt 

since tt is Lipschitz continuous. In order to deal with the next integral, we observe that 

\F\iPa,i3) - F'{ipfi)\ < a.e. in Q, 


due to the mean value theorem and fl2.49p . Hence, applying the Holder and Sobolev 
inequalities, and owing to estimate fl2.60p for ipa^g and the regularity 02.391) of ipp, we can 
infer that 


/ {F'{p>a,l3) - F'{pi3))fl < C -F |(p/3p -F ijlg 

JQt Jo 

<C [ (||7'«,/3|l6 + W^aWl + 1) llT'lh yhds 


: ds 


10 


< c 


>0 


H\\p\\vds<6 [ \\fi\\yds + cs [ |v3p. 

Jo JQt 


Next, we take advantage of 04.Sp for dtjj.a,p and deduce that 

-a [ + fdp) ds < 6 [ {\\ip\\y + \\n\\v) ds + csa^ [ 

Jo Jo Jo 

— (llv^lly + llhlly) ds + csa. 

Jo 

Finally, the last term in 04.61) can be simplihed and estimated in the following way: 


P - P)iP + (3p-(j)+ 


' Qt 


<p {ap - p(p + {/3 + l)ap) <6 |/ip cs 


' Qt 


' Qt 


' Qt 


+ kP). 


By combining 04.61) and these inequalities, choosing 5 small enough and applying the 
Gronwall lemma, we obtain 02.50p . 
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